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ABSTRACT 

We present a complete 12 CO J = 1 — ■> map of the prominent molecular ridge in the Large 
Magellanic Cloud (LMC) obtained with the 22-m ATNF Mopra Telescope. The region stretches 
southward by ~ 2° (or 1.7 kpc) from 30 Doradus, the most vigorous star-forming region in the Local 
Group. The location of this molecular ridge is unique insofar as it allows us to study the properties of 
molecular gas as a function of the ambient radiation field in a low-metallicity environment. We find 
that the physical properties of CO-emitting clumps within the molecular ridge do not vary with the 
strength of the far-ultraviolet (FUV) radiation field. Since the peak CO brightness of the clumps shows 
no correlation with the radiation field strength, the observed constant value for CO-to-H2 conversion 
factor along the ridge seems to require an increase in the kinetic temperature of the molecular gas 
that is offset by a decrease in the angular filling factor of the CO emission. We find that the difference 
between the CO-to-H2 conversion factor in the molecular ridge and the outer Milky Way is smaller than 
has been reported by previous studies of the CO emission: applying the same cloud identification and 
analysis methods to our CO observat ions of the LMC molecular ridge and CO data from the outer 
Galaxy survey by I Dame et al.l (|2001|) . we find that the average CO-to-H 2 conversion factor in the 
molecular ridge is Xco — (3-9 ± 2.5) x 10 20 cm~ 2 (K km s -1 ) -1 , approximately twice the value that 
we determine for the outer Galaxy clouds. The mass spectrum and the scaling relations between the 
properties of the CO clumps in the molecular ridge are similar, but not identical, to those that have 
been established for Galactic molecular clouds. 

Subject headings: Magellanic Clouds — ISM: molecules — ISM: structure - galaxies: ISM 



1. INTRODUCTION 

The study of the interstellar medium (ISM) in dwarf 
galaxies provides important insights into the star forma- 
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tion history of the Universe. In massive spiral galax- 
ies the differential rotation causes shear forces that in- 
fluence the morphology and evolution of the progenitor 
molecular gas and the young stellar component, making 
the study of the their relationships difficult. In general, 
dwarf galaxies show little evidence of shear so we can 
study the relation between star formation and its prereq- 
uisites at their very locations. They may also be consid- 
ered, in a hierarchical galaxy formation scenario, as local 
templates of the earliest galaxies. They are thought to 
be survivors of the merging processes that lead to the for- 
mation of larger galaxies, such as the Milky Way. Thus, 
the detailed study of the low-metallicity molecular gas 
in dwarf galaxies may help us to understand the physi- 
cal conditions under which star formation took place in 
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the main population of galaxies in the early Universe. 
This understanding is crucial in order to interpret obser- 
vations of the ISM in galaxies at high redshift. 

In dwarf galaxies, stars form out of low- 
met ahicity J _^ow_djasMoj : gas ratio dense molecular 
gas (jLisenfeld fc Ferraral fl998h . Far-ultraviolet (FUV; 
6eV < hv < 13.6 eV) photons from newly formed 
massive stars illuminate the progenitor molecular 
clou ds producing photon-dom inated regions (PDRs; 
e.g. iHollenbac h fe Tielendll999l . and references therein). 
As the FUV field is mainly attenuated by dust ab- 
sorption, the penetration of FUV photons is larger in 
low-metallicity molecular gas than in its solar metal- 
licity counterpart. This produces an enhancement 
of the photodissociation rate of CO (increasing the 
relative abundance of C and C + ) while H2 remain s 
almost unaffected (e.g. Ivan Dishoeck fc Blackl Il98l . 
The relations between properties of the molecular gas 
(e.g. mass, CO luminosity, CO-to-H 2 conversion factor) 
with metallicity and the strength of the FUV radiation 
field have been investiga t ed from bot h theo r etical (e.; 
Malonev fc Blackl 119881: lElmegreenl Il989t ISakamotc 



19961) and empiric al 
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Rubio et all 119911: IWilsonl Il995t lArimoto et alj 11996. 
Israeli 119971 |2000D perspectives. The CO-to-H 2 conver- 
sion factor (X co = iV(H 2 )//co; an empirical factor 
that relates the intensity of CO [ico] to the column 
density of H2 LV(H2)]) has been found to increase for 
clouds located in galaxies wi th low met allicities and 
strong FUV radiation fields ([Israel 1993). However, 
Xqo h as also been show n to decrease with the FUV 
field by IWeiB et al.l (|2001[ ). This apparent contradiction 
arises because Xqo depends on both the metallicity and 
the FUV radiation field, but mutual variation of these 
parameters in different interstellar environments makes 
it difficult to determine how they each influence the 
value of Xqo 7 ■ Additional complications are introduced 
when low-metallicity molecular clouds are observed on 
different spatial scales: in particular, Xco has been 
found to vary wit h the resolution of the observations 
(|Rubio et alJ 119931 : IVerter fc Hodgd Il995f ). At higher 
angular resolution, the measured properties of molecular 
clouds in low metal licity galaxies tend to approach those 
of Galactic clouds | Johansson et a.1.1 Il998h iBolatto et al.1 
120031: iLerov et all 120061 ). Ideally, we would like to 
conduct high resolution studies of molecular clouds in 
systems where one property of the interstellar envi- 
ronment varies while other environmental conditions 
remain constant. 

The Large Magellanic Cloud (LMC) is an excellent ex- 
tragalacti c laboratory for such a study. Its proximity 
(~50 kpc: lFeast|[l999l a value which we will use through- 
out the paper) and nearly fa ce-on orientation (i = 35°; 
Ivan der Marel fc Cionil [200D allows us to study individ- 
ual clouds at high spatial resolution, while strong spatial 
variations in the LMC's far-ultraviolet (FUV) radiation 

7 Note tha t Xco a ^ so depends on the gas density and kinetic 
temperature l|Malonev fc B lack 1988). However, the 12 CO J = 
1 — * line is likely to be thermalized at the typical densities of 
molecular clouds and therefore its emission is not sensitive to the 
gas density. The kinetic temperature is related to the FUV field 
because the heating in the region where CO is emitted is dominated 
by photoelectric heat ing, which depends on the FUV field (e.g. 
Bakes & Ticlcns 1993). 



field provide an opportunity to study the influence of 
the local FUV field strength on the prope rties of low- 
meta llicity molecular gas (30 — 50% solar; Wes terlundl 
[1991 . 

The most striking feature in the LMC is the 30 Doradus 
H 11 complex (hereafter 30 Dor) , which is the most mas- 
sive star-forming region in the Local Group. Pre- 
vious survey s of molecular g as in the LMC (e.g. 
Cohen et all Il988t iFukui et al.l fl999l Uohansson et al.l 



1998; iKutner et al.l ll997) have revealed a prominent com- 



plex of molecular clouds extending south from 30 Dor 
for nearly 2 kpc, which we refer to as "the 30 Dor 
molecular ridge" . Studies of individual clouds in the 
30 Dor molecular ridge have found gas densities between 
~10 4 — 10 5 c m" 3 and kinetic temperatures between 
~10— 80K (e.g. [Johansson et al.lll99 8: Mina midani et Id] 
l2007l:lPineda et al.ll2008l) ~ 

The young stellar cluster R136 located at the cen- 
ter of 30 Dor produces an extre me FUV radiatio n 
field. As seen in radio continuum (jDickel et al.l [2 005^ . 
Ha dGaustad et al.ll200l . far-infrared (FIR; |Miwering| 
119891 ) and ultra-violet (e.g. ISmith et al.lll987l ) maps, the 
strength of the field decreases with increasing distance 
from R 136. The 30 Dor molecular ridge thus presents an 
ideal interstellar environment in which to investigate the 
properties of low-metallicity molecular gas as a continu- 
ous function of the ambient FUV radiation field. 

Another aspect of interest is the global structure of the 
low-metallicity interstellar medium (ISM) of the LMC. 
Cloud structure is often characterized by the clump-mass 
spectrum, which describes how the total cl oud mass is 
distrib uted among different-size structures. IFukui et a.1.1 
(|2001h presented a mass spectrum for GMCs in the 
LMC, based on NANTEN 12 CO data with spatial res- 
olution of ^40 pc. These authors derived a power-law 
index a = 1.9 for the LMC mass spectrum, which is 
similar to that found in the Galaxy (1.4 to 1. 9 ; e.g . 
Elmegreen fc Falgaronel 119961: iHeithausen et all 119981 : 



Kramer et alj|1998t [Williams et al.lll994D . Whether this 



relation holds at even smaller scales, where clumps have 
low masses and are therefore more susceptible to the ef- 
fect of enhanced CO photo-dissociation in low-metallicity 
environments, is still unclear. This is mainly due to the 
lack of wide-field, high resolution maps of 12 CO emission 
from molecular clouds in low metallicity environments. 

In this paper, we present a complete, fully sampled, 
12 CO J = 1 — > map of the 30 Dor molecular ridge in 
the LMC obtained with the ATNF Mopra 22-m Tele- 
scope. Parts of this region have been mapped with the 
SEST telescope at a similar angular resolution and in 
the s ame transition ([Johansson et al.lll998t IKutner et al.l 
Il997f h Cloud properties such as sizes, CO luminosities, 
Xco factor, etc. ha ve been studied i n the 30 Dor and 
N159/160 regions bv I Johansson et aTl (|1998f) but n ot in 
the southern molecular clouds (jKutner et a.1.1 Il99~7h . al- 
though the gravitational stabil ity of clumps in t his re - 
gion have been discussed by llndebetouw et ail (|2008h . 
We combine our 12 CO data with a map of the strength 
of the FUV radiation field, derived from HIRES/IRAS 
60 /im and 100 /im images, in order to study the de- 
pendence of the properties of low-metallicity molecular 
clouds on the radiation field. We also present a study 
of the clump-mass spectrum for the molecular clouds in 
the 30 Dor molecular ridge. This has not been attempted 
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in the LMC at the resolution of our observations. The 
main difference between our work and that done with the 
SEST telescope is in the analysis of the data. Addition- 
ally, the wider coverage of our map allows us to identify 
clouds that hav e not previously been mapped at this res- 
olution (see also lOtt et aU 2008). The paper is organized 
as follows. The observations are summarized in Sect. O 
The analysis of the data is described in Sect. [3] and the 
results are discussed in Sect. fJJ A summary is given in 
Sect. [5] 

2. OBSERVATIONS 

The 12 CO J = 1 -> (115 GHz) observations were 
made using the 22-m ATNF Mopra telescope between 
June and September 2005. The 12 CO J = 1 -> in- 
tegrated intensity map of the 30 Dor molecular ridge is 
shown in Fig. [lj The o bservational resu lts of this map 
are presented elsewhere (|Ott et al.ll2008f ). 

A total of 120 fields, 5'x5' in size, were mapped using 
the on-the-fiy mapping mode, encompassing the entire 
30 Dor molecular ridge. Typical values of the single side- 
band (SSB) system temperature T* ys of about 600 K were 
obtained under clear atmospheric conditions. Prior to 
each map, the pointing accuracy was checked using the 
SiO maser R Dor, with corrections typically below 5". 

The digital autocorrelator was configured to output 
1024 channels across a 64 MHz (163.8 km s -1 ) band in 
two polarizations, corresponding to a velocity resolution 
of 0.16 km s" 1 . We used the AIPS++ task LiveData 
to calibrate the spectra and to remove a first-order fit, 
and another AIPS++ task Gridzilla to grid the data 
into spectral line cubes. LiveData is the processing soft- 
ware originally designed for data taken with the Parkes 
multibeam receiver and is used to apply system temper- 
ature calibration, bandpass calibration, heliocentric cor- 
rection, spectral smoothing, and to write out the data in 
sdfits format. Gridzilla is a re-gridding software pack- 
age that forms three dimensional (RA-Dec- velocity) data 
cubes from bandpass-calibrated sdfits files (usually pro- 
duced by LiveData). During the gridding procedure, the 
data were convolved with a Gaussian smoothing kernel 
with a FWHM half of the beam-size (33" at 115 GHz), 
producing an effective resolution of 45" for the data cube. 
The spectra were smoothed to an angular resolution of 
60" (corresponding to a linear scale of ~15 pc) and to a 
velocity resolution of 0.96 km s^ 1 , in order to increase the 
signal-to-noise ratio. The Mopra telescope has an "in- 
ner" error beam in the range of 40 — 80", which is similar 
to the size of emission. We the refore applied an e xtended 
beam efficiency of ?7 m b = 0.55 (|Ladd et al.ll2005l ) to con- 
vert from T£ to extended beam brightness (T xrj ). The 
typical r.m.s. noise of the channel maps is 0.16 K. 

In order to compare the properties of the molecular 
clouds located in the 30 Dor molecular ridge with Galac- 
tic clouds of approximately solar metallicity, we utilize a 
map of the 2nd Quadrant observed with the CfA 1.2m 
telescope and presented by Dame etaUlpOOl . The 2nd 
Quadrant map covers the Galactic plane from I — 100° 
to 150° and b = -2° to 2°, with a noise level of about 
0.1 K, at a resolution of 450". We use this dataset (kindly 

8 The Mopra radio telescope is part of the Australia Telescope 
which is funded by the Commonwealth of Australia for operation 
as a National Facility managed by CSIRO. 
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Fig. 1. — Contour map of the integrated intensity of the 12 CO 
J = 1 — » line observed toward the 30 Dor molecular ridge. 
This contour map is overlaid onto an Spitzer 24 /an emission map 
(Mcixncr ct al. 2006). To improve the contrast, we divide the con- 
tour map into white and black regions. The contour levels corre- 
spond to the 5% to 85% of the peak intensity (53 K km s — 1 ), in 
steps of 10%. 
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provided by Thomas Dame) because it samples Galactic 
clouds at physical scales that are similar to those of our 
30 Dor molecular ridge map. A further advantage of this 
data is that blending of different clouds along each line 
of sight is negligible for outer Galaxy clouds. While the 
properties of molecular clouds in the outer Galaxy are 
known to be dif ferent to those in the inner Galaxy (e.g. 
iMav et al.|[l997ft . another reason to compare our molec- 
ular ridge data to the outer Galaxy cloud sample is that 
molecular clouds in the inner Milky Way are subject to 
conditions that characterize the central regions of spiral 
galaxies (e.g. shear forces, tidal forces, strong magnetic 
fields, cosmic rays), and which may not be common in 
the LMC. 

3. ANALYSIS AND RESULTS 
3.1. Clump Decomposition 

In order to study the properties of the molecular gas in 
the 30 Dor molecular ridge and compare them to prop- 
erties of clouds in the outer Galaxy, we identify individ- 
ual clumps in the 30 Dor molecular ridge and 2nd Quad- 
rant datasets using the software package G AUSSCLUMPS 
(jStutzki k. Gustenlll990HKramer et al.lll998| ). By apply- 
ing the same decomposition algorithm to both datasets, 
we ensure that differences between the properties of 
molecular clouds in the two samples are not due to the 
decomposition algorithm. GAUSSCLUMPS decomposes the 
observed three-dimensional (position and velocity) inten- 
sity distribution into emission from individual clumps 
by fitting Gaussians to individual peaks and subtracting 
them in an iterative process. The iterative decomposition 
process is terminated when the peak intensities of a new 
generation of clumps fall below 5 times the r.m.s. noise 
of the original datacube. The decomposition process 
reports clump parameters such as their position, Local 
Standard of Rest (LSR) as velocity, clump orientation, 
spatial FWHMs along the spatial axes (Ax, Ay), bright- 
ness temperature (T m h), and FWHM linewidths (Av). 
iKramer et al.l (|1998f) described in detail the stability of 
the results to variations of input parameters used in the 
decomposition of several molecular clouds. For the de- 
composition of our observations, we adopt the standard 
set of stiffness parameters used in their analysis. 

As GAUSSCLUMPS approximates irregularly-shaped fea- 
tures by Gaussian clumps, the algorithm can report the 
identification of small, low-mass clumps that are the ac- 
tually the residue of previously subtracted larger struc- 
tures. To reduce the impact of this effect on our anal- 
ysis, we only accept Gaussian clumps with fitted sizes 
that are larger than 110% of the spatial resolution (i.e. 
the deconvolved size of the clump is at least 40% of the 
instrumental resolution along both spatial dimensions). 
The FWHM linewidth of a clump is limited to be larger 
than 2.44km s _1 and 1.63km s _1 for the 30 Dor molec- 
ular ridge and 2nd Quadrant datasets respectively, i.e. 
we reject clumps that occupy a single pixel in velocity 
space. At the observed spatial scales, we do not ex- 
pect that this criterion excludes a significant number of 
ge nuine clumps wit h smal l Av. The catalog presented 
by iJohansson et al.l (|1998t ) has only one clump with a 
FWHM linewidth less than 2.44 km s _1 , even though 
their velocity resolution is ~0.1 kin s . Moreover, only 
four clumps with linewidths narrower than 1.63 km s _1 
are identified in the 2nd Quadrant (a similar number is 



seen in the catalog presented bv lMav et al.lll997t h 

In order to determine the kinematic distances to the 
identified clumps in the 2nd Quadrant dataset, we as- 
sume a flat Galactic rotation curve with rotation velocity 
©lsr = 220 kms -1 and radius Rq = 8.5 kpc. We dis- 
carded all (194) clumps with kinematic distances lower 
than 2 kpc, as the estimation of kinematic distances is 
not accurate enough for such local emission. 

We identify a total of 67 and 118 clumps in the 30 Dor 
molecular ridge and 2nd Quadrant datasets, respectively. 
The average (deconvolved) radius of the identified clumps 
in the 30 Dor molecular ridge and 2nd Quadrant is 
~10pc, with a dispersion of ~ 35%. 

The derived properties of the identified clumps in the 
30 Dor molecular ridge are listed in Table 1. The CO 
luminosity of a clump is estimated using 
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We derive the virial mass, M v - lr , using 
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where R is the deconvolved clump radius and Aw' is the 
deconvolved FWHM of the observed line profile. The fac- 
tor k depends on the a ssumed density profile of a clump 
(jMacLaren et al.l [T988) . In the following, we adopt a 
value oik = 210, corresponding to a spherical cloud with 
uniform density. Different values of k (i.e. different den- 
sity profiles) only change the scaling but not the slope 
of any relation involving M v j r . The deconvolved clump 
radius is defined as 
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pc 



Ax' 
pc 



1/2 



pc 
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where Ax' and Ay' are the projected deconvolved 
FWHM angular sizes of the clump. We choose to use de- 
convolved values in the calculation of the virial mass be- 
cause it is an intrinsic property of the cloud. Some pub- 
lished studies do not deconvolve the cloud parameters, 
and instead quote cloud property measurements that are 
convolved with the instrumental resolution. This leads 
to an overestimate of the virial mass, in particular if the 
cloud size and FWHM linewidth are of the same order as 
the instrumental resolution. For example, the values of 
Xqo obtained in our survey (see Section l3~4f are ~ 40% 
larger if we use undeconvolved values for the clump ra- 
dius and FWHM linewidth; for the 2nd Quadrant data, 
using undeconvolved values increases the derived Xqo 
values by ~ 30%. The different methods used to calcu- 
late the virial mass and other cloud parameters compli- 
cates a quantitative comparison between our results and 
those of previous studies. 

We estimate the value of Aco for clumps in the 30 Dor 
molecular ridge assuming that they are in virial equi- 
librium. In this case, Aco is calculated dividing the 
chimp's virial mass - which we assume is an accurate 
measure of the total mass of H2 present - by its CO 
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FlG. 2. — Cumulative histograms of the FWHM linewidth and 
CO luminosity resulting fro m our decomposit i on of t he 30 Dor ridge 
and the data published by I Johansson et al.l Jl998ft . The FWHM 
linewidth is not deconvolved by the instrumental resolution. For all 
curves, the horizontal error bars reflect the uncertainties in the ob- 
served properties and the vertical error bars represent y/W counting 
errors. 

luminosity. These quantities are the result of multi- 
plying the peak H2 column density and CO intensity 
(Ljo) of a (Gaussian) clump by its area. Therefore, 
X co = N(R 2 )/I co = M vir /L CO - Note that L co does 
not need to be calculated from deconvolved values of R 
and Av since Lco is independent of the resolution of the 
observations. 

The effect of adopting different values of k in eq. 
is to reduce Xco by 10% and 66%, for density profiles 
of the form n(r) oc r _1 and n(r) oc r~ 2 , respectively. 
However, the relative values of Xco m the 30 Dor ridge 
and the 2nd Quadrant do not depend on the adopted 



density profile, as their virial masses were calculated with 
the same formula. 

The errors in the clump properties are estimated by 
considering the r.m.s. noise in the positions where 
clumps were identified and by assuming that the clump 
parameters obtained by GAUSSCLUMPS (Aa;, Ay, Av) have 
errors of 30%. We use error propagation to obtain uncer- 
tainties for Xco j Mvin and Lco- The virial masses and 
CO luminosities have typically errors of 50%, whereas 
the error of Xco is typically 70%. 

We test our decomposition by comparing our measure- 
ments of the properties of clumps in th e 30 Dor molecu- 
lar rid ge to those listed in Tab le 1 of I Johansson et all 
(1998). Uohansson et alj (|1998f ) derive cloud parame- 
ters by fitting circular Gaussians to each intensity peak 
in their integrated intensity maps, so an exact corre- 
spondence between the results of the two studies is not 
to be expected. In Figure [2l we show cumulative his- 
tograms of the FWHM linewidths 9 and the CO luminosi- 
ties of clumps identified by us and by I Johansson et al.l 
(1998). We exclude clumps in the 30 Dor molecular 
ridge dataset that are located outside the 30 Dor an d 
N159/160 regions observed by Uohansson et al.1 (|1998h . 
Furthermore, w e only consider clumps identified by 
iJohansson et alj (|1998l ) that have similar CO luminosi- 
ties to the clumps in the 30 Dor molecular ridge, as their 
observations are more sensitive than ours. Note that 
we choose Aw and Lco because they are independent 
of angula r resolution which is d ifferent in the 30 Dor 
ridge and IJohansson et al.l (|1998l ) datasets. For our de- 
composition, the mean (istandard deviation) Lco is 
(6.9±7.9) x 10 3 K km s" 1 pc 2 and Av is (4.9±2.1) kms" 1 . 
In case of IJohansson et all (|1998l ). the mean (istandard 
deviation) L C o is (6.6±8.0)xl0 3 Kkms _1 pc 2 and Av 
is (5.3±1.9 kms" 1 ). Kolmogorov-Smirnov (KS) tests 
between the distributions of Av and Lco for the two 
datasets indicate that the probabilities (P) that the cloud 
samples are drawn from the same parent distribution 
are 0.2, and 0.1, respectively (We consider P< 0.05 to 
be statistically significant and P< 0.1 to be marginally 
significant). The Lco and Av measurements for the 
IJohansson et a l. ( 1998) and 30 Dor ridge datasets cover a 
similar range of values, ev en though differences in decom- 
position methods used bv IJohansson et ail (|1998h and us 
almost certainly affect the exact distribution of the clump 
properties. 

In the following, we compare the results from our de- 
composition of the 30 Dor molecular ridge dataset with 
those of the decomposition of the 2nd Quadrant data. 
As noted above, this ensures that any difference (or sim- 
ilarity) between these two datasets is not an effect of the 
adopted decomposition method. 

3.2. Estimation of the Far-ultraviolet Radiation Field 

An important feature of the 30 Dor molecular ridge 
is that the low-metallicity molecular gas is exposed to 
very different strengths of the FUV radiation field. This 
can be seen in maps of different tracers of the FUV 
radiation field, e.g. at radio continuum (jDickel et al.l 

9 Note that we use undeconvolved values of the FWHM 
linewidth for th i s com parison, in order to be cons istent with 
IJohansson et af] l|1998T l . As explained in Section 13.41 we use de- 
convolved values of the clump parameters in our analysis. 
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TABLE 1 

Clump properties in the 30 Dor molecular ridge 



Clump 


q(J2000) 


<5(J2000) 




Av' a 


R a 




log(L co ) b 


log(M vir ) 6 




x c 


# 


[h m s] 


[°'"] 


[kms- 1 ] 


[kms- 1 ] 


[PC] 


[K] 


[K km s- 1 ] 


[M ] 


10 20 [cm- 2 /(K km s" 1 )] 


[xo] 


1 


05:36:14.455 


-69:01:48.77 


272.1 


7.6±2.3 


7.6±1.6 


1.3±0.2 


.3.7 


5.0 


8.1±5.8 


5 


2 


05:38:30.816 


-69:02:04.52 


249.6 


4.6±1.4 


8.7±1.8 


1.6±0.1 


3.7 


1.6 


3.6±2.5 


119 


3 


05:35:56.511 


-69:02:17.78 


259.4 


5.3±1.6 


10.8±2.3 


2.3±0.1 
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a Deconvolved by the instrumental resolution. 

The CO luminosities and virial masses have typically errors of 50%. 
C Strength of the FUV radiation field in units of the Drainc i 1978) field. They have typical errors of 20%, as this is the accuracy of the HIR.AS/IR.AS images used to derive this 
quantity (see text). 
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Fig. 3. — Strength of the FUV radiation field as a fun ction of 
distance to 30 Dor. The values of \ are m units of the IDraine! 
1978 field. Each data-point represents the value of \ at the central 
position of each identified clump. The main and secondary peaks 
corresponds to clumps close to the 30 Dor and N159 regions. The 
lines represents the geometrical dilution of the FUV field in case 
30 Dor and N159 are considered the only sources of FUV radiation 
(see text). 



20051) . Ha (jGaustad et al.ll2001l), far-infrared fSchwerind 



1989), and ultra-violet (e.g. ISmith et all Il987f ) wave- 
lengths. For all tracers, the intensity peaks in regions 
that are closely associated with star formation, and de- 
clines in locations that are more distant from regions of 
active star formation. Assuming that metallicit y gradi- 
ents are negligible in the 30 Dor molecular ridge (|Dufourl 
[1981 . we can therefore study clump properties as a func- 
tion of the strength of the FUV radiation field. To de- 
rive a map of the strength of the FUV radiation field 
in the 30 D or molecular ridg e , we follo w the method 
presen ted in iMookeriea et al I (|2006l ) and iKramer et al.1 
(2008]). A map of the FIR continuum intensity (/fir) is 
constructed using HIRES/IRAS 100 /im and 60 fim im- 
ages obtained from the IPAC data center. The resolution 
of these images has been enhanced to about 110"and 
72", respectively usin g a maximum correlation method 
l|Aumann et al.lll990l ). The 60 /mi image is then con- 
volved to 110", in order to match the resolution of the 
100 /im image. To estimate /fir between 42.5 pjn. and 
122. 5 Atm, we use the relat ionship (|Helou et al.lll988i see 
also iNakagawa et all I1998D : 



/fir 



erg cm 2 s 1 sr 



3.25 x 10~ u 



l 



Jysr" 1 



1.26 x 10~ 



/lOOAin 



Jysr 1 



(4) 



The dust grains in molecular clouds absorb most of the 
incident FUV photons, which are later re-radiated in the 
FIR continuum. Assuming that the contribution from 
other heating mechanisms such as cosmic rays can be 
neglected, the strength of the FUV field can be derived 
from /fir according to 



_\_ 

Xo 



4tt 



/fir 



1.71 (2 • 1.6 x lO- 3 ) erg cm- 2 s- 1 sr" 1 ' 



(5) 



Here, x lS expressed in units o f the iDrainel (|1978l ) 
field (xo), which is related to the [Habing) (|1968l ) field 



(1.6x10 erg s cm ) by a factor of 1.71 when av- 
eraged over the full FUV spectrum. It is likely that the 
contribution of FUV photons wit h hv < 6 eV increases 
the d ust heating by a factor of ~ 2 (|Tielens fc Hollenbachl 
1985), and that the bolometric d ust continuum in tensity 
is a factor of ^2 larger than /fir (|Dale et al.ll2001l ): these 
two corrections approximately cancel out. Note that in 
the case where molecular clouds do not fill the beam, 
this method for estimating the strength of the FUV field 
provides only a lower limit. 

Typical angular sizes of clumps are between ^66" to 
200", and as the resolution of the x ma P is about 110", 
its value represents an average over an entire clump. The 
calibration of the x map appears to be consistent with 
previous estimates of the FUV radiation field in the vicin- 
ity of the 30 Dor molecular ridge. Near the N159W cloud 
(Clump #22) for example, we obtain x— 173, which is 
in good agreement with an independent determination 
from UV, radio continuu m, and IRAS 60 /x m emission 
( X = 120 - 350) made by llsrael et all (|1996f) . Our esti- 
mate for x m N159W is also consistent with the results 
of PDR models for this m olecular cloud (e.g. iPak et aTl 
Il998t iPineda et all l2008h . Although the x=173 may 
be considered low for a mass ive star-forming region (in 
the Orion Bar, for exa mple, iWalmslev et al. | |2000l and 
lYoung Owl et all 120001 find y ~ 10000). llsrael et all 
(|1996| ) argue that the low metallicity and low dust-to- 
gas ratio of this region increase the FUV photon free 
path lengths, resulting in a greater geometric dilution of 
the FUV radiation field. We stress, however, that for 
this analysis we are interested in the dependence of the 
CO clumps parameters on relative variations of the FUV 
radiation field. 

In Fig. [31 we show x as a function of the distance 
to 30 Dor. Each point corresponds to the value of x 
in the central position of each molecular clump identi- 
fied in the 30 Dor molecular ridge. The value of x for 
each clump is also listed in Table 1. The strength of the 
FUV field peaks in clouds near the 30 Dor and N159 Hn 
regions, and then decreases for more distant positions. 
The thick line represents the expected distribution of x 
if 30 Dor and N159 were the only sources of FUV radi- 
ation, where we have used the value of x a t the clump 
closest to the center of each H ll region to represent the 
radiation source, and considering only geometrical dilu- 
tion of the FUV field. The individual contribution from 
each of these regions is shown by a dashed line. Fig. [3] 
shows that the radiation from 30 Dor and N159 is not 
sufficient to explain the observed distribution of \, and 
that additional OB stars along the 30 Dor ridge must 
therefore contribute to the local FUV radiation field. 

3.3. Clump properties as a function of the FUV 
radiation field 

Figure S] shows Lqo, M V - 1T , and Xqo plotted as a func- 
tion of the FUV radiation field. The horizontal lines 
represent the mean values of the data given in Table 2. 
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Fig. 4. — The CO-to-H2 conversion factor (Xqo), the virial mass 
(M v j r ), and the CO luminosity (Lqo) as functions of the strength 
of the FUV radiation field. The horizontal lines represent the mean 
values of the data (see Table 2). 

As we can see, these quantities do not show any corre- 
lation over three orders of magnitude in x- To interpret 
the plots of -M v i r and Leo with we show the bright- 
ness temperature (T m b), clump FWHM linewidth (Au'), 
and clump radius (R) as a function of \ m Figure 
Again, the horizontal lines represent the mean values of 
the data given in Table 2. These basic clump param- 
eters also appear to be independent of the strength of 
the FUV radiation field. In Section 01 we discuss possi- 
ble reasons why the properties of the CO clumps in the 
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Fig. 5. — The brightness temperature (T mb ), clump FWHM 
linewidth (At/), and clump radius (Ft) as functions of the strength 
of the FUV radiation field. The horizontal lines represent the mean 
values of the data (see Table 2) 

30 Dor molecular ridge are insensitive to the strength of 
the FUV field. 

3.4. Comparison with clouds in the outer Galaxy 

We compare the properties of clumps in the 30 Dor 
Ridge with those of the 2nd Quadrant Galactic data. 
The mean values and standard deviation of these prop- 
erties are listed in Table 2 and their cumulative his- 
tograms are shown in Figures and [7] Although mea- 
sures of their dispersion and central tendency are com- 
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parable, Kolmogorov-Smirnov (KS) tests between the 
distributions of Av', R, and T m b for the two datasets 
indicate that the probabilities (P) that the cloud sam- 
ples are drawn from the same parent distribution are 0, 
0.07, and 0, respectively. These low probabilities sug- 
gest that there are differences between the distribution 
of CO clump properties in both regions. KS tests for the 
distributions of £co, M v i T , and Xqo also result in low 
probabilities (P=0.3, 0, and 0), mainly due to the dis- 
similarity between the distribution of Av', T m b, and R 
in the 30 Dor ridge and the 2nd Quadrant. In Figures [6] 
we can see that clumps in the 30 Dor Ridge have lower 
brightness temperatures compared with clumps in the 
2nd Quadrant. This might be an indication that the CO 
beam filling factor is lower in the 30 Dor Ridge, which can 
be a result of its reduced metallicity. The 30 Dor Ridge 
clumps have larger FWHM lincwidths: 50% of clouds in 
the 2nd Quadrant have FWHM linewidths less than 4 
km s _1 compared to only 30% in the 30 Dor Ridge. This 
is in ag r eemen t with the predictions from the theory of 
iMcKeel (| 19891) (see Section 3.4.1). Finally, the 30 Dor 
Ridge also shows a deficit of clumps with R < 8 pc com- 
pared to the 2nd Quadrant. This might be the effect of 
CO photodissociation that destroy or prevent the forma- 
tion of small clumps. The differences in Av' , T m b, and R 
result in lower CO luminosities and larger virial masses 
in the 30 Dor Ridge. We caution that in case of Av' and 
R the observed variations are comparable to the errors 
in the determinations of these quantities. 

The mean value of the Xqo conversion factor in the 
30 Dor molecular ridge is 3.9±2.5xl0 20 cm~ 2 (K km 
s -1 ) -1 ). This is almost twice the average value for 
X co in the 2nd Quadrant ([2.0±1.0]xl0 20 (cm" 2 K km 
s -1 ) -1 ), although we note that the dispersion of X C o 
values for the clumps in both datasets is quite large. Ta- 
ble 2 and Figure [7] show that the average virial mass of 
clumps in the molecular ridge is ~ 45% larger than the 
average virial mass of the 2nd Quadrant clouds, but that 
the average CO luminosity is ~ 45% lower for the molec- 
ular ridge clumps. This differences account for the factor 
of ^2 larger value of Xqo m the 30 Dor molecular ridge. 
The average value of Xco m the 30 Dor molecular ridge 
is consis tent with th e 4.3 x lO 20 cm~ 2 (K km s -1 ) -1 de- 
rived by llsrael et al.l (|2003j) using SEST observations of 
several molecular cloud associated with Hn regions in 
the LMC (not including clouds in the 30 Dor Ridge). 

3.4.1. Empirical Relations 

Several studies of molecular clouds in our Galaxy have 
shown that there are empirical scaling relations between 
the basic observable physical properties of molecular 
clouds (e.g. [La rson 198lt ' (Mverdll983t iDame et al.lll986t 
ISolomon et al.|[l987l ). These relations, often referred to 
as "Larson's Laws" , link the FWHM linewidth and size 
of a cloud, its CO luminosity and FWHM linewidth, its 
virial mass and CO luminosity, and its CO luminosity 
and radius ( note that only two of t hese relations are in- 
dependent). ISolomon et alj (|1987l ). for example, found 
that clouds in the inner Milky Way follow power-laws 
of the form Av oc R 5 , L co oc Aw 5 , M vir cx L co l , 
and Lco °c i? 2 ' 5 . Power-laws with similar slopes have 
already been observed for LMC molecular clouds, with 
some evidence that clouds in the LMC have lower CO 
luminosities compared to Galactic clouds with a similar 
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Fig. 6. — Cumulative histograms of the brightness temperature 
(T m b), clump FWHM linewidth (Av'), and clump radius (R) de- 
rived in the 30 Dor Ridge and 2nd Quadrant. For all curves, the 
horizontal error bars reflect the uncertainties in the observed prop- 
erties and the vertical error bars represent %/~N counting errors. 



FHW M linewidth (e.g. ICohen et alj Il988t iGarav et all 
Il993t Uohansson et alj|1998t iMizuno et al.ll2001a| ). This 
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Fig. 7. — Cumulative histograms of the virial mass (M v i r ), the 
CO luminosity (Lqq), and the CO-to-ti2 conversion factor (X^q) 
derived in the 30 Dor Ridge and 2nd Quadrant. For all curves, 
the horizontal error bars reflect the uncertainties in the observed 
properties and the vertical error bars represent V~N counting errors. 

has been previously interpreted as evidence that CO is 
significantly photo-dissociated in the LMC. 



TABLE 2 

Mean Values of Derived CO Clump Parameters 



Parameter Mean STDV 



Units 



30 Dor Ridge 



^CO 

M vir 

£co 
R 

Av' 



3.9 
5.2 
7.6 
10 
4.6 
1.7 



2.5 
4.7 
8.3 

3 

1.8 
0.9 



10 a) cm- 2 (K km s" 1 )" 1 
10 4 Mq 
10 3 K km s _1 pc 2 
pc 
km s — 1 
K 



2nd Quadrant 



^CO 
M vir 

Leo 
R 

Av' 

Tmb 



2.0 
3.6 
11 
10 
3.8 
3.0 



1.0 
3.9 

18 
3.8 
1.4 
1.6 



10 2u cm- 2 (K km s" 1 )" 1 
10 4 M 
10 3 K km s" 1 pc 2 
pc 
km s~ 1 
K 



Figured] shows the observed relations between size, ve- 
locity dispersion, CO luminosity and virial mass for the 
clumps identified in the 30 Dor molecular ridge and in the 
2nd Quadrant. The virial mass-CO luminosity relation 
is shown in Figure^. The orthogonal BES fit 10 results 



M vir 



_ 10 (0.8±0.4) r (1.0±0.1) 



J CO 



for the 30 Dor molecular ridge, and 



M v 



_ in (1.3±0.2) r (0.8±0.04) 
— IU 1j co 



(G) 



(7) 



for the 2nd Quadrant. The relation is slightly steeper 
for clumps in the 30 Dor molecular ridge compared to 
the molecular clouds in the 2nd Quadrant, such that the 
larger (M > 5 x 10 4 M Q ) clumps in the 30 Dor molec- 
ular ridge appear to be underluminous by a factor of 
~2.4, relative to a 2nd Quadrant cloud with a similar 
virial mass. In Section 13.51 we will see that the num- 
ber of massive clumps is larger in regions with stronger 
FUV radiation field (x > 10) than in regions with weaker 
FUV radiation field (x < 10). This suggests that effects 
of CO photodissociation (enhanced by the FUV radia- 
tion field) might be more significant in high-mass rather 
than in low-mass clumps. This might explain the appar- 
ent underluminosity of clumps with large virial masses 
in the 30 Dor Ridge. We note that the best-fitting virial 
mass-CO luminosity relation for clumps in the 30 Dor 
molecular ridge is linear, which implies that the value of 
-^co (°c M v i r / Lqo) is independent of the clump mass. 
In the 2nd Quadrant, by contrast, the value of Xqo de- 
creases with increasing mass. The variation in Xco is 
quite small, however, corresponding to a factor of 0.3 de- 
crease in Xco over roughly two orders of magnitude in 
mass. 

Figured shows the FWHM linewidth-CO luminosity 
relation. The fits are 

L GO = 10(2-2±O.3) At/ (2.3 ± O.5) (g) 

for the 30 Dor ridge, and 

L CO = 10(l-7±°^) Au '(3.8±0.4) (g) 

for the 2nd Quadrant. The 30 Dor ridge shows a shal- 
lower slope compared with that for outer Galaxy clouds. 

10 We fit power-laws using the orthogonal bi-variate e rror and in- 
trinsic scatter method (BES; Akritas &i Bcrshadv 1996). Available 
at www . astro .wise . edu/~mab/archive/stats/stats .html. 
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Fig. 8. — The virial mass-CO luminosity, CO luminosity-FWHM linewidth, CO luminosity-clump size, and FWHM linewidth-sizc 
relations for the identified clumps in the 30 Dor molecular ridge. For comparison, we also include the same relations for the clumps 
identified in the 2nd Quadrant. The thick and thin lines correspond to fits to the relations for the 30 Dor ridge and 2nd Quadrant, 
respectively. For clarity, we do not show error bars. The typical errors in virial mass, CO luminosity, FWHM linewith, and size are 50%, 
50%, 30%, and 20% respectively. 

for the 2nd Quadrant. Both 30 Dor Ridge and the 2nd 
Quadrant show a similar slope. This slope is close to 
Lqo R 2 ' 5 , which is expected for clouds in virial equilib- 
rium that follow M vir oc Lqq and Av' oc R 05 . The best- 
fitting CO luminosity-size relations that we determine 
imply that the mass surface density of the clumps in- 
creases slowly with increasing clump mass, although the 
average mass surface density for the large clumps in both 
datasets is in good agreement with the canonical value for 
molecular clouds in the inner Milky Way (170 M Q pc~ 2 , 



Clumps in the 30 Dor ridge with large FWHM lincwidths 
(Aw' > 3 km s _1 ) have CO luminosities that are a factor 
of 1.4 lower than those of clumps in the 2nd Quadrant 
with similar Av'. 

The CO luminosity-size relation is shown in FigurelHfc. 
The fits are 

L co = 1 (1-3±0-4)J 2 (2.4±0.4) ( 1Q ) 

for the 30 Dor ridge, and 

L CO = l (l-4±0.2) i? (2.5±0.2) (n) 
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assum ing the average Xco values that were estimated 
above iSolomon et all I1987L but see IHever et alJ (|2009l ) 
for evidence that this should be revised downwards to 
- 100 Mq pc- 2 ). 

Unfortunately, we are not able to fit a power-law to 
the FWHM linewidth-size relation shown in Figure[8]i, 
as the dynamic range in Av' and R is insufficient. Nev- 
ertheless, the close agreement between the mean Xco 
derived from our a nalysis of the CO em ission and the 
value determined bv lDobashi et al.1 ((2008) from dust ex- 



tinction measurements suggests that the clumps in the 
30 Dor molecular ridge are indeed close to virial equi- 
librium. As the clump properties do not show any sys- 
tematic variation along the ridge (Figure 2]), we conclude 
that the degree of virialization is uniform for all clumps 
(see also llndebetouw et al.1 120081 ). Although, we see in 
FigureEJi that the clumps in the 30 Dor ridge and the 
2nd Quadrant occupy the same region in the parameter 
space, Table 2 and Figure [HI suggests that clumps in the 
30 Dor Ridge have slightly larger Av' relative to the 2nd 
Quadrant. This is in apparent agreement with the theory 
of phot o ioniza tion-regulated star formation proposed by 
iMcKeel (fl989), which predicts that molecular clouds in 
low metallicity galaxies will have larger velocity disper- 
sions than clouds of a similar size with solar metallicity. 
Nonetheless, wc note that this prediction assumes that 
the filling fraction of the CO emission, the densities of 
the CO-emitting clumps and the clump-to-cloud extinc- 
tion ratios are similar for molecular clouds in different 
galax ies, which may not be the case (e.g. IHever et al.1 
l200l . 

In conclusion, we find small differences between the 
scaling relations exhibited by the clump properties in 
the outer Galaxy and in the 30 Dor molecular ridge. This 
suggests that the basic physical properties of the molecu- 
lar clumps are not especially sensitive to variations in the 
metallicity of the interstellar gas, or that the difference 
in metallicity between the 30 Dor molecular ridge and 
clouds in the outer Galaxy is too small to produce signif- 
icant disparities between the clump properties. A similar 
result h as been obtain e d acro ss a wider range of metallic- 
ities bv iBolatto et all (|2008l ) (see also lLerov et al.ll2006l 
and lSheth et al.l 120081) . 

3.5. Clump-Mass Spectrum 

A general property of molecular clouds is that they 
show self-similar structure over a wide range of scales. 
This hierarchical, clumpy structure may be quantified 
through the clump-mass spectrum of molecular clouds, 
which describes how the molecular gas is distributed 
across low- and high-mass substructures. Velocity- 
chan nel maps of molecu l ar clouds in the Milky Way 
e.g. IWilliams et allll994lElmegreen fc Falgaronelll996l; 



(e.g. [W illiams et al .lliiJiMt Ib lmcgrccn & t algaroneyiyytt 
Stutzki et al.lll988HKramer et al.lll998l ; lHeithausen et all 



1998; IHever fc Terebevl I1998D and in external galaxies 
( Sta nimiroyic et apl2000HStanimirovic fc Lazarianf2001l ; 
iFukui et al.ll200H iKim et al.ll2007ft can be decomposed 
into clumps that exhibit a power-law mass spectrum of 
the form dN/dM oc M a , where the exponent, a, is typ- 
ically observed to lie between —1.4 to —1.9. This value 
of a, which indicates that the molecular gas is prefer- 
entially distributed within massive substructures, holds 
over several orders of magnitude of clump mass (e.g. 
iHeithausen et al1ll998l ). and is independent of the st.ar- 




M "i: ■: (M e ) 

Fig. 9. — Cumulative clump-mass spectra for the (upper panel) 
30 Dor ridge and (lower panel) the 2nd Quadrant. The lines over 
the clump spectra represents the best truncated power-law fit over 
a range of masses above the completeness limit (M > 1.76xlO 4 M0; 
vertical line). 



forming nature of the clouds (e.g. [Williams et"al"1ll994f ). 

In this Section, we study whether the metallicity and 
FUV radiation field modify the distribution of clump 
masses in the ISM of the LMC. For this analysis, we cal- 
culate the clump- mass spectrum of the 30 Dor molecular 
ridge. We also calculate the mass spectrum of the clumps 
identified in the 2nd Quadrant, in order to compare the 
mass spectrum of clumps in low-metallicity and solar- 
metallicity gas, avoiding systematic effects produced by 
the selection of a particular decomposition method. 

In the literature, the clump-mass spectrum based on 
12 CO data is often calculated using clump masses de- 
rived from the clump luminosity by applying a uniform 
Xqo conversion factor. Since the scatter in Xco is large 
(Section l3.4[) . in this study we prefer to use virial masses 
for the calculation of the mass spectrum. The total virial 
mass of the 30 Dor molecular ridge is about 3.5 x 10 6 M . 

We express the clump-mass spectrum in its cumulative 
form rather than in its differential form. This is moti- 
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TABLE 3 

Parameters of Mass Spectra for clumps in the 30 Dor 
ridge and 2nd quadrant 



Dataset 


a 


N 


Af„/(10 5 M Q ) 


30 Dor ridge 


-2.0 ± 0.4 


5.9 ± 5.9 


2.5 ± 0.9 


2nd Quadrant 


-2.0 ± 0.3 


5.8 ± 6.2 


2.4 ± 1.1 


30 Dor ridge (x > 10) 


-1.8 ± 0.3 


1.3 ± 2.7 


5.5 ± 1.6 


30 Dor ridge (x < 10) 


-1.2 ± 0.3 


97 ± 10 


1.2 ± 0.3 



vated by the relatively low number of identified clumps 
(N = 67) in the 30 Dor ridge. For such a low number 
of clumps, the slope of the differential clump-mass spec- 
trum is sensitive to the binning process (which is not 
required by the cu mulative mass spectrum). Following 
IRosolowskvl (f2005t ). we fit a truncated power law of the 
form 



N(M' > M) = N 



M 
Mo 



- 1 



(12) 



to our cumulative clump mass spectrum, where Mq is 
the maximum mass of the distribution and N is the 
number of clumps with masses larger than the mass 
where the distribution shows significant deviation from 
a power-law, given by 2 1 ' a+1 M$. A truncation at high 
masses in the cumulative mass spectrum can either be a 
real effect (i.e. due to a physical mechanism that regu- 
lates cloud growth) or th e result of finite sampling (e.g. 
iMcKee k Williams! [l997l) . 

We use the erro r-in-variables procedure described by 
IRosolowskvl (|2005l ) to fit a power-law to cumulative mass 
spectra. This approach has the advantage that it ac- 
counts for both the (AN) 1 / 2 statistical error in the distri- 
butions and the errors in the measurements of the clump 
masses. 

The introduction of a cutoff criterion that rejects noise 
features and clouds with low CO emission in the decom- 
position implies that some low-mass clumps are missing 
from the clump-mass spectrum. This effect produces a 
turnover in the mass spectra at low masses. We assume 
that the incompleteness of our sample begins to be signif- 
icant when the clump mass is lower than 5 times the min- 
imum mass that can be selected in our sample based on 
the criteria described in Section 13.11 for the 30 Dor ridge 
(3.5 x 10 3 M Q ). We therefore fit a truncated power law 
to the cumulative clump mass distribution using clumps 
with masses greater than our adopted completeness limit 
of ~ 2 x 10 4 M . For consistency, we use this com- 
pleteness limit for both the 30 Dor molecular ridge and 
the 2nd Quadrant datasets. 

Fig. [5] shows the cumulative clump-mass spectra of 
the 30 Dor molecular ridge and the 2nd Quadrant. The 
results of the fits are given in Table 3. 

The exponents of the power-law fits to the clump- 
mass spectra of the 30 Dor ridge and in the 2nd Quad- 
rant are similar within the errors. This suggests that 
the reduced metal abundance in the LMC does not sig- 
nificantly affect the clump-mass spectrum of molecu- 
lar clouds across the observed range of clump masses. 
The fits are also consistent w i th the mass spectra that 
were derived by IRosolowskvl ((2005) for clouds in the 
outer Milky Way (a = -2.29±0.08, 7V =4.5±3.5, and 




id 3 io 4 10 5 id 6 

Mass >M,y 

Fig. 10. — Cumulative clump-mass spectra for the 30 Dor ridge 
for clumps with \ > 10 (upper panel) and with \ < 10 (lower 
panel). The lines over the clump spectra represents the best trun- 
cated power-law fit over a range of masses above the completeness 
limit (M > 2 X 10 4 M o ; vertical line). 



M O /1O 5 M =2.9±1.O; data from lBrunt et al.ll2003D and 
in the entire LMC (a = -1.71 ±0.19, iVo=10±6.5, and 
Af o /1O 5 M =23±4.6; data from lMizuno et al.ll2001ri ) us- 
ing the same fitting procedure and virial estimates for 
the molecular c loud masses. Not e that the truncated 
power-l aw fit by [Rosolowskyj (|2005T ) to the iMizuno et al.l 
(|2001bf ) LMC cloud catalog was based in virial mass mea- 
surements that do not account for beam deconvolution. 
If we fit a truncated power law to our molecular ridge 
data using undeconvolved measurements for the clump 
properties, we obtain a slope for the mass spectrum of 
a = —2.2 ± 0.3, which is still close to their results. Some 
of the scatter in the derived slopes is probably due to 
the different decomposition algorithms used to identify 
cloud structures in the various studies, although we note 
that the cloud mass spectrum constructed from the LMC 
cloud catalogue includes clouds across the LMC, not just 
the molecular ridge region. It also extends to larger cloud 
masses and imposes a higher completeness limit than we 
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have adopted for our analysis. Previous works that find 
flatter slopes for the clum p mass spectrum in the outer 
Galaxy (e.g. a = -1.45 iMav et al.l I1997D are not di- 
rectly comparable to our results, as the derived slope is 
strongly dependent on how the mass spectrum is con- 
structed (jRosolows kv 2005). For this reason, we prefer 
to compare our results with independent analyzes that 
have employed the same representation of the mass spec- 
trum that we have used here. 

We verified that the use of virial masses instead of 
masses derived from CO luminosities and a constant 
value for the Xco conversion factor does not have a sig- 
nificant impact on the clump mass spectra of the 30 Dor 
molecular ridge and 2nd Quadrant. For clumps in the 
molecular ridge, assuming Xco=3-9xl0 20 cm -2 (K km 
s -1 )^ 1 (see Section(2H]), the best-fitting truncated power 
law to the clump mass spectrum yields a = —2.2 ± 0.3, 
N =1.5±2.8, and M /10 5 =2.4±0.7. For the 2nd Quad- 
rant, the best fit gives a = —1.9 ± 0.2, N =5±8, and 
M o /10 5 =4.8±2. 7, assuming X co =2. Ox 10 20 cm" 2 (K km 
s -1 ) -1 . For both datasets, the fits are consistent with 
those determined for the cumulative mass spectra con- 
st ructed from yirial e stimates of the clump masses. 

iFukui et al.l (|2001f ) obtain a slope for the molecular 
cloud mass spectrum in the entire LMC of —1.9, covering 
a range of cloud masses from ~ 8 x 10 4 to 3 x 10 6 M Q . 
Our analysis shows that the power-law relation extends 
to even lower masses (~ 2 x 10 4 M Q ), suggesting that 
the ISM in the LMC might show the self-similarity over 
a wide range of sc ales, as observed for mo lecular clouds 
in the Milky Way (|Heithausen et al.|[l998l ). 

In order to study the influence of the strength of the 
FUV radiation field on the clump mass spectrum, we 
divide the clumps in the 30 Dor molecular ridge into low- 
FUV (x < 10) and high-FUV (x > 10) subsamples. The 
resulting clump-mass spectra are shown in Figure fTOl and 
the results of the truncated power-law fits to the spectra 
are listed in Table 3. The derived slope for the x < 10 
subsample (-1.2±0.3) is dominated by low mass clumps, 
and the mass spectrum steepens considerably for larger 
clump masses. If we fit the spectrum with a simple 
power-law (i.e. without truncation) of the form 

N(M> > M) = ^— J , (13) 

where Mo and a represent the same quantities used 
in Equation (| 1 2[) . we obtain a much steeper slope of 
a = — 2.7 ± 0.2. Figure fTOl shows that there is a deficit 
of high-mass clumps with \ < 10. This might be be- 
cause massive stars form preferentially in massive rather 
than in low-mass clumps. Since the strength of the local 
FUV field depends on embedded massive stars, the aver- 
age FUV field will be larger for massive clumps than for 
low-mass clumps, leading to the apparent deficit of large 
clumps with % < 10. 

4. DISCUSSION 

As outlined in the previous sections, the properties of 
molecular clouds in the 30 Dor molecular ridge, as traced 
by their CO emission, are not significantly affected by the 
varying intensity of the FUV radiation field interacting 
with the reduced- metallicity gas of the LMC. The requi- 
site physical conditions for CO emission may provide a 
possible explanation for this lack of dependence. As the 



photodissociation of CO is enhanced in low-metallicity 
systems, CO emission is only present in regions with 
large H2 column densities. Thus, CO observations might 
only trace regions that are wel l-shielded from the FUV 
radiation, as suggested by e.g. Uohansson et al.1 (|1998l) . 
In this case, the 30 Dor molecular ridge might contain a 
substantial fraction of H2 coexisting with other carbon 
species such as C and C + , and this "CO-dark" fraction 
of the molecula r gas may be large r than for clouds in 
the Milky Way (iGrenier et al.l l2005l 10 to 50%). Indeed, 
iPineda et al.l (|2008f ) recently found that the column den- 
sity of CO toward LMC N159W is very much reduced 
compared to the column densities of C and C+. This 
might be a reason why studies of molecular clouds in 
low-metallicity galaxies based on dust emission (which 
should trace the entire molecular content of clouds once 
the contribution from H I is subtracted) usually mea- 
sure larger m olecular masses than studies based on CO 
emiss i on (e.g. [I srael 19971 iRubio et al.ll2004t iLerov et al.1 
l2007t iBot et al.l I2007D . However, the measurement of 
molecular gas masses from dust emission usually assumes 
that the dust and gas are well mixed, and further relies 
on estimates of the dust-to-gas ratio, dust temperature, 
and optical properties of dust grains. Using near infrared 
extinction, a method that is independent of the cloud's 
dynamical sta t e and of cloud decomposition method, 
IDobashi et~ail (|2008l ) obtained an average Xco conver- 
sion factor of 2.26xl0 20 cm~ 2 (K km s -1 ) -1 in the en- 
tire LMC which is a factor of ~1.7 lower than our mean 
value of the 30 Dor ridge 11 . Note that the IDobashi et all 
(2008) value refers to an H2 column density associated 
w ith CO-emi t ting g as . 

Elmcgrecn ( 1989) modelled the relation between Xco 
and the strength of the FUV field as X co oc x 3/S /T, 
where T is the 12 CO J = 1 — > brightness tempera- 
ture. We can test this prediction using our da ta (see Fig- 
ure Hp . Note that the relation presented by lElmegreenl 
(1989) also depends weakly on the clump mass and on the 
ambient pressure, although we ignore their contribution 
here. We corrected the observed 12 CO J — 1 — > peak 
bright ness temperature assuming a beam filling factor of 
1/6 (|Pineda et al.l |2008). Figure ITU shows that our data 
is i nconsistent (y 2 ~ 254) with the theoretical prediction 
by lElmegreenl (|1989l ). Our result that Xco is not corre- 
lated with x for the molecular clumps in our survey also 
appears to c ontradict the Xco oc X relation obtained by 
llsraell (|1997f) . However, the range of FUV radiation field 
estimates in that paper (which are derived from infrared 
observations and expressed in terms of the radiative en- 
ergy per nucleon) is quite limited, whereas our estimate 
of x extends over about three orders of mag nitude. We 
also note that the cloud mass estimates from llsrael ()1997f ) 
are based on dust observations. If Xco is independent 
of x for CO observations and Xco X f° r dust obser- 
vations, this would imply that the amount of "CO-dark" 
molecular gas increases with the strength of the FUV 
radiation field. However, thi s would seem to be in con- 
tradiction with the results of IDobashi et al.l ((2008), who 

The value of Xco derived from dust extinction depends on 
the adopted extinction to total column densit y ratio. When us- 
ing the values adopted by Imara & Blitz (2007), their Xco factor 
is 3.55xl0 20 cm -2 (K km s — in closer agreement with our 
determination. 
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Fig. 11. — Relation between Xco, T and the strength of the 
FUV radiation field for the clumps in the 30 Dor ridge dataset 
(solid black boxe s). The straight line is the relation predicted by 
lElmeereenl lTl989h . 

found no systematic variation of the extinction-based 
X co as a function of distance to 30Doradus (i.e. from 
quiescent regions to regions exposed to strong FUV ra- 
diation fields). 

The observed relation bet ween Xco a nd y can be rec- 
onciled with the theory of Elmcgrcen (1989) if the ki- 
netic temperature scales with the FUV radiation field. 
Since the CO brightness temperature does not show sig- 
nificant variations with x (Figure [5]), this would require 
that the beam filling decreases with x throughout the 
30 Dor Ridge. It is plausible that the kinetic is larger 
in regions with higher \, since photoelectric heating is 
likely to be the dominant heating mechanism. For ex- 
ample, if clouds with % = 1 have T^ n = 10 K, kinetic 
temperatures of 32 K and 66 K are required for N159W 
(x = 173) and 30Dor-10 (x = 407), respectively, to com- 
pensate any increase of Xco with x- These kinetic tem- 
peratures are consistent with estimations in the litera- 
ture Ce.fi. [Johansson et alJll998t iMinamidani et alJl2507t 
iPineda et al J 12008). A systematic decrease of the beam 
filling with x might be the result of an enhanced rate of 
CO photodissociation which reduces the size of CO cores 
within the beam. However, other mechanism might also 
play a role e.g. cloud disruption by stellar winds from 
newly formed stars, etc. 

The average Aco conversion factor found in our sur- 
vey of the 30 Dor molecular ridge is a factor of ~ 2 
larger than that found in outer Galaxy clouds, a differ- 
ence that is smaller compared with previous studies of 
clouds in the LMC at coarser angular resolutions. As we 
can see in Fig[TJ the molecular ridge is composed of a few 
cloud complexes with linear scales between 100 — 300 pc. 
These cloud complexes show substructure characterized 
by compact clumps connected by low-surface brightness 
gas on scales between 10 - 30 pc. Since GAUSSCLUMPS 
decomposes the CO emission and identifies clumps by 
fitting local CO intensity peaks, this substructure is the 
one characterized in the present paper. Therefore, the 
average Xco estimated in our survey refers to scales be- 



tween 10 — 30 pc. The large cloud complexes do not nec- 
essarily have the same Vco as the individual clumps. 
Although increases in the clump size are tracked by 
increasing estimates for the virial mass, the CO lumi- 
nosity will not show an equivalent increase if a signifi- 
cant quantity of low-surface brightness inter-clump gas 
is present. In this case, the appropriate Vco conver- 
sion factor will be higher for larger scale molecul ar gas 
structures, as shown by e.g. iBolatto et alJ ([2003) . Ob- 
servations at coarser angular resolution trace these large- 
scale complexes rather than small clumps, which helps 
to explain the larger values of Vco derived by previous 
surveys in the LMC. An example of how Vco increases 
with the resolution of the observations is shown in Fig- 
ure 1121 where we plot the average Vco obtained when 
GAUSSCLUMPS is applied to versions of the 30 Dor molecu- 
lar ridge datacube that have been smoothed to different 
angular resolutions. We include the average Aco for 
2.6'(40pc) and 8.8' (127 p c), corresponding to the res- 
olutio n of the NANTEN (|Fukui et all [2008) and CfA 
1.2m (Co hen et al.lll98"ct ) telescopes. At the resolution 
of the NANTEN survey, we obtain an average value of 
Vco ~5.1xl0 20 cm -2 (K km s -1 )" 1 , which is somewha t 
lower than the value obtained by iFukui et al.l (2008), 
Vco = 7 x 10 20 cm~ 2 (K km s" 1 ) -1 . Our result is in 
better agreement with the determination by iBlitz et al.l 
(2007) for the same dataset, Aco = 5.4 x 10 °cm 
CK km s" 1 )" 1 : the discrepan cy between the IBlitz et al.1 
|2007t ) and lFukui et~aT] ([20081 ) may arise from a different 
choice of input paramet ers for the cloud decompos ition 
algorithm (CL0UDPR0PS; R osolowskv fc Lerovll200ll . At 
the resolution of the ICohen et al.1 (|1988f ) observations, 
we obtain A C o — 6.9x 10 20 cm -2 (K km s " 1 )" 1 , a fac- 
tor of ~1.6 lower than the value quoted by ICohen et al.1 
(1988) (V co - 11 x 10 20 cm~ 2 (K km s" 1 )" 1 , assuming 
1.8xl0 20 cm- 2 (K km s" 1 )" 1 for Milky Way clouds). In 
summary, Figure [12] suggests that scale of the observed 
structures must be considered when the total molecu- 
lar mass of galaxies is calculated and that comparisons 
between extragalactic and Galactic molecular cloud pop- 
ulations should adopt the same decomposition and cloud 
identification methods. 

5. SUMMARY AND CONCLUSIONS 

We have presented a large-scale 12 CO J = 1 — ► 
survey of the 30 Dor molecular ridge in the Large Mag- 
ellanic Cloud observed with the Mopra telescope. From 
the analysis of these observations we conclude the follow- 
ing: 

1. We estimate the strength of the far- ultraviolet ra- 
diation field for the entire 30 Dor molecular ridge, 
in order to study the dependence of CO cloud prop- 
erties on this parameter. We find that CO clump 
properties are insensitive to changes in the strength 
of the FUV radiation field over almost 3 orders of 
magnitude in x- The observed relation between 
CO-to-H2 conversion factor and the strength of the 
FUV field can be reconciled with the expected in- 
crease predicted by theory if the kinetic temper- 
ature also increases in regions with strong FUV 
fields. Since the CO brightness temperatures do 
not show any significant variation along the 30 Dor 
molecular ridge, this would require that the angu- 



16 



Pineda, Ott, Klein, Wong, Muller & Hughes 



RKoLlion (Arc m in) 




SO 100 
R«oLrlion (Paraec) 

Fig. 12. — Average Xco obtained for the 30 Dor ridge map 
for different angular resolutions. The error bars correspond to the 
standard error of the mean. 

lar filling factor of the CO emission shows a corre- 
sponding decrease with \. 

2. Applying the same cloud decomposition method to 
clouds in the 30 Dor molecular ridge and in the 
outer Galaxy clouds in the 2nd Quadrant, the av- 
erage CO-to-H 2 conversion factor obtained in our 
survey is Xqo — 2Imw- This discrepancy is 
smaller than has been reported by previous studies 
of the CO emission in the LMC at lower angular 
resolution. 

3. Studying the scaling relations between the physical 
properties of the observed clumps, we find that the 
molecular clouds in the 30 Dor ridge exhibit prop- 
erties that are similar, but not identical, to those 
in the outer Galaxy. In particular, we find some 
evidence that clumps in the 30 Dor molecular ridge 



have larger velocity dispersions and lower bright- 
ness temperatures than their outer Galaxy coun- 
terparts. Additonally, we find a deficit of small 
clumps in the 30 Dor molecular ridge compared to 
the outer galaxy. Further studies of the emission 
from atomic hydrogen and other carbon species in 
the interclump gas within LMC molecular clouds 
would be useful for determining the physical origin 
of these slight variations. The degree of virializa- 
tion of the clumps does not appear to show any 
significant variation along the molecular ridge. 

4. An analysis of the clump mass spectra yields a 
power-law index of a = —2.0 ± 0.4 for the entire 
30 Dor molecular ridge. The clump mass spectrum 
is similar to the mass spectra obtained for Galac- 
tic clouds and for the entire LMC cloud catalog. 
Comparing different sub-regions with different av- 
erage strengths of the FUV field, we find that the 
clump-mass spectrum is shallower in regions that 
are exposed to strong FUV fields. We propose that 
this is caused by the preferential formation of mas- 
sive stars within large molecular clumps, rather 
than an effect of CO photo-dissociation in the low- 
metallicity molecular gas. 
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